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Abstract 

The AdS/CFT correspondence relates Wilson loops in M = 4 SYM 
theory to minimal area surfaces in AdS§ space. In this paper we con- 
sider the case of Euclidean flat Wilson loops which are related to 
minimal area surfaces in Euclidean AdS% space. Using known mathe- 
matical results for such minimal area surfaces we describe an infinite 
parameter family of analytic solutions for closed Wilson loops. The 
solutions are given in terms of Riemann theta functions and the valid- 
ity of the equations of motion is proven based on the trisecant identity. 
The world-sheet has the topology of a disk and the renormalized area 
is written as a finite, one-dimensional contour integral over the world- 
sheet boundary. An example is discussed in detail with plots of the 
corresponding surfaces. Further, for each Wilson loops we explicitly 
construct a one parameter family of deformations that preserve the 
area. The parameter is the so called spectral parameter. Finally, for 
genus three we find a map between these Wilson loops and closed 
curves inside the Riemann surface. 
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1 Introduction 



One of the first results of the AdS/CFT correspondence [JJ was the com- 
putation of Wilson loops and from there the quark anti-quark potential as 
done by Maldacena, Rey and Yee [2J. Although much work was devoted to 
the computation of Wilson loops only few explicit examples are known of 
minimal area surfaces in AdS 5 space. In the case of closed Euclidean Wilson 
loops (with constant scalar) the most studied one is the circular Wilson loop 
[3] which is dual to a half-sphere. The only other one we are aware of is 
the two intersecting arcs (lens shaped) [4j. For infinite Wilson loops, parallel 
lines [2] and the cusp are known [5]. In the case of multiple contours as for 
example two concentric circles, interesting results were found using integra- 
bility [6]. In the language we use here they correspond to elliptic functions 
which appear for genus one g = 1. In the case of Minkowski signature AdS 
space and in particular light-like lines more is known starting with the light- 
like cusp [7] and culminating with a large recent activity [5] in relation to 
scattering amplitudes following [91 [TO]. 

In this paper we point out that in the case of flat Euclidean Wilson loops 
which are dual to minimal area surfaces in Euclidean AdS 3, much can be 
done by using known results from the mathematical literature [UJ. In fact, 
an infinite parameter family of solution is known in terms of Riemann theta 
functions. This type of construction using theta functions is described in 
detail in [12] and was already used in the case of strings moving in t x S 3 
by Dorey and Vicedo [13] and in the case of an Euclidean world-sheet inside 
AdS 3 space by Sakai and Satoh in [H]. Here we consider Euclidean Wilson 
loops inside Euclidean AdS 3 and rederive the original results by perhaps 
more pedestrian methods based on the trisecant identity for theta functions. 
In this way theta functions are thought as special functions whose properties 
fit well with the equations of motion of the string in AdS 3 space much in 
the same way as trigonometric functions fit the harmonic oscillator equation. 
Each theta function and therefore each Wilson loop is associated with an 
auxiliary Riemann surface of given genus g. A relatively simple formula is 
derived for the area and an example for genus g = 3 is worked out in detail. 
Perhaps our main contribution is to find closed Wilson loops and to derive 
a formula for the renormalized area that follows the AdS / CFT prescription. 
The calculations are done at the classical level, it should be interesting to 
extend them for example to one-loop as can be done in the case of the circular 
Wilson loop [T5] . 
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This paper is organized as follows. We start by writing the equations of 
motion and use the Pohlmeyer reduction procedure to simplify the equations 
and arrive at the cosh-Gordon equation plus a set of linear equations. In the 
following section we review the properties of the theta functions and show 
how they can be used to solve the equations of motion and compute the 
regularized area. Finally we construct a particular example of genus three 
where we show that there are closed Wilson loops that can be described by 
this method. We plot the corresponding surfaces and compute the areas. Be- 
sides, we also describe the mapping of the Wilson loop into a curve embedded 
inside the Riemann surface. In the last section we give our conclusions. 



2 Equations of motion 

In this section we write the equations of motion and simplify them using 
the Pohlmeyer reduction [16]. In the context of Minkowski space-time this 
procedure was used by Jevicki and Jin [17J to find new spiky string [18] 
solutions and by Alday and Maldacena [10J to compute certain light-like 
Wilson loops. In the case of Euclidean AdS 3 that we are interested in here, 
we can use embedding coordinates X /J=0 ...3 parameterizing a space _R 3,1 and 
subjected to the constraint 

Xl-Xl-Xl-Xl = l, (2.1) 

with an obvious 50(3,1) = SX(2,C) global invariance. The space has an 
S 2 boundary at infinity. Other useful coordinates are Poincare coordinates 
(X, Y, Z) given by: 

X + zY= Xl + lX2 , Z = ^^. (2.2) 
Xo — X3 Xo — X3 

The boundary is now an R 2 space and located at Z = 0. A string is pa- 
rameterized by world-sheet coordinates a a = (a, r) or equivalently complex 
coordinates z = o + ir, z = a — ir. The action in conformal gauge is given 
by 

S = l -J (dX ll dX' M -A{X lt X ,t -l)) dadr (2.3) 
= l -J ^{d a Xd a X + d a Yd a Y + d a Zd a Z) dadr, (2.4) 
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where A is a Lagrange multiplier and the \i indices are raised and lowered 
with the R 3 ' 1 metric. An Euclidean classical string is given by functions 
Xp(z, z) obeying the equations of motion: 

dBX^ = AX M , (2.5) 

where A, the Lagrange multiplier is given by 

A = -dX^dX". (2.6) 

These equations should be supplemented by the Virasoro constraints which 
read 

dX^dX" = = BX^BX^. (2.7) 

Later on we will be interested in finding the solutions in Poincare coordi- 
nates (X, Y, Z) but for the moment it is convenient to study the problem in 
embedding coordinates X^. We can rewrite the equations using the matrix 

X =U°^ 2 x~-t )=W. <«> 

where a % denote the Pauli matrices. Notice also that Poincare coordinates 
are simply given by 

Z=^, X + zY = ^. (2.9) 

-^22 ^22 

The matrix X satisfies 

X f =X, detX = l, <9<9X = AX, det(dX) = = det(<9X) , (2.10) 



as follows from the definition of X, the constraint (2.1), the equations of 
motion (2.5) and the Virasoro constraints (2.7). We can solve the constraint 
XJ = X by writing 

X = AA t , detA = l, Ag5L(2,C). (2.11) 

The equations of motion have a global SL(2, C) = SO (3, 1) symmetry under 
which 

X -»■ UXU\ A -)■ UA, U e SL(2, C). (2.12) 
In the new variable there is an SU(2) gauge symmetry 

A^AU, IA (z, z) G SU (2) , (2.13) 



4 



since this leaves X invariant. We can define the current 

J = A- 1 dA, J = A~ 1 dA, (2.14) 

which is invariant under the global symmetry and, under the local symmetry 
transform as 

J ^UUU + U ] dU, J ->■ U ] JU + U ] dU. (2.15) 

From the definition of J, J, the property that detA = 1, the equations of 
motion and the constraints we find: 

dJ-dJ+[J,J] = 0, (2.16) 
TrJ = TrJ = , (2.17) 

<9(J + J f ) + \[J-J\J + J f ] = 0, (2.18) 

det(J+J t ) = 0, (2.19) 

and the corresponding equations found by hermitian conjugations of the ones 
given. In the third equation we used that for SU{2) currents we have for 
example: 

JJ= 1 -[J,J] + 1 -Tt(JJ) , (2.20) 

and similarly for the other products. The trace part gives the Lagrange 
multiplier A as: 

A = ^Tr((j+Jt)(jt + J)) , (2.21) 

which does not provide an equation but is useful later to determine the 
world-sheet metric. From the form of the equations it seems convenient to 
define 

A=\(J + f), B=\(J-P). (2.22) 
The equations read now 



Tr^l = TiB = , (2.23) 

det.4. = , (2.24) 

dA+[B,A] = 0, (2.25) 

dB + dB ] = [B\B] + [A\A]. (2.26) 



5 



The SU(2) gauge symmetry acts on these currents as 

A^U ] AU, B^U ] BU + U ] dU, U(z, z) e SU{2). (2.27) 

In a sense, B plays the role of a gauge field. Since TiA = we can write in 
terms of Pauli matrices a- 5 : 

A = (Uf + iAf) a j , (2.28) 

with A 1 ' 2 two real three-dimensional vectors. The property det .4 = implies 
that they are orthogonal and of the same length. The SU(2) gauge symmetry 
acts on them as three-dimensional rotation so we can take A^ to be along 
the x axis, and A^ 1 along the y axis. In that way we can choose the gauge 
such that 

A = I e «M)( ai + lG2 ) = e ^) a+ ; (2.29) 

where a(z,z) is a real function and a + = +icr 2 ) = ^ q q J • Equation 
dA + [B, A] = uniquely implies that 

B = -h)aa z + b 2 (z, z)a + , (2.30) 

for some function 6 2 (^,^)- Finally the equation for B implies that b 2 = 
f(z)e~ a for an arbitrary holomorphic function f(z). It also implies that 

dda = e 2a + ffe~ 2a . (2.31) 

So, up to a gauge transformation the most general solution is given by 

A = eV+ , (2.32) 
B = - l -daa z + f(z)e- a a + , (2.33) 



with a satisfying eq.(2.31). Finally the function f(z) can be locally elimi- 
nated by changing to world-sheet coordinates w such that dw = ^ffdz. If 
we further redefine a. — > a + \ \n(ff) and do a gauge transformation with 

U = e~2 az (where i<p = \ f/f) then the result is equivalent to setting / = 1. 
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From the equations of motion for A and B, namely eqns.(2.23)-(2.26) it 



can be seen that they are invariant under multiplying A by a constant of 
modulus one that we call A (|A| = 1). We get: 



A 
B 



AeV 
1 



daa 7 



dda = 2 cosh(2a) 



(2.34) 
(2.35) 
(2.36) 



where we set / = 1 by the reasons indicated before. The constant A can be 
eliminated by a gauge transformation and a redefinition of a but we keep it 
for later convenience. It is called the spectral parameter and should not be 
confused with the coupling constant in the dual gauge theory. 



Having computed A and B we can use eq.(2.22) to reconstruct J and J 
obtaining: 



J 



J 



— \da 



— e 



(2.37) 



Finally we should use eq.(2.14) to compute A. Summarizing we need first to 
solve the equation: 

dda = 2 cosh 2a , (2.38) 



then plug a into the definitions for J, J, namely eq.(2.37), and solve for 



<9A 
BA 



AJ , 
AJ. 



(2.39) 
(2.40) 



Finally, the string solution is determined as X = AA\ The equation for a is 
non-linear but the ones for A are linear since J, J are known once a is known. 
This is the main idea of the Pohlmeyer reduction [16] which we rederive here 
as it applies to our particular problem. Similar considerations in the context 
of string theory are well-known, for example see [17], [19], [10], j!4j . 

Notice that, since TrJ = TrJ = we automatically find that detA is 
constant independent of z, z. However we need det A = 1 so we just need to 
normalize A dividing by an appropriate constant. Furthermore it is conve- 
nient to write 

4>i i>2 



(2.41) 
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where the vectors ip — (^1,^2) an d ip = (^1,^2) are linearly independent 
and satisfy 

dip = ipj, dip = ^J , (2.42) 

and the same for if>. They have to be linearly independent so the determinant 
(■01^2 — 4>24>i) is non vanishing (but is constant as discussed before). Even 
with these conditions there is a certain ambiguity in choosing ip, ip but those 
boil down to SL(2, C) = 50(3, 1) transformations of X. 



3 Solutions 



As shown in [TT] solutions to eqns. d2~38|2.39|2.40| ) can be found using theta 



functions. In this section we rederive the results of [TT] using a perhaps more 
pedestrian way. The method we use is to consider the theta functions as 
special functions whose derivatives are such that they are good candidates to 



solve eqns.( 2^38|2.39|2.40 ) in much the same way as trigonometric functions 
are good candidates to solve the harmonic oscillator. The equations are 
solved by simple substitution and adjustment of the parameters. To start we 
need to review some properties of the theta functions. 

3.1 Riemann theta functions and their properties. 

There is a vast literature on Riemann theta functions [20J. In this section we 
review the minimal knowledge necessary to find solutions to the equations. 
We follow the notation of [2T] which also gives a good introduction to Rie- 
mann surfaces. Notice also that in the next section we develop an example 
in detail which can be read in parallel with this section. Consider a compact 
Riemann surface of genus g with fundamental cycles dj, bi {i = 1 . . . g) and 
intersections 

<2j o a,j = = bi o bj, a, o bj = 5ij , (3.43) 

which means that only intersects 6j. The Riemann surface is taken to be 
a hyperelliptic one defined by the function: 



MA) 



29 

A;Q(A-A,). (3.44) 
\ 3=1 

The square root has cuts with branching points at 0, 00, Xj but is well defined 
in a double cover of the complex plane. This double cover is the Riemann 
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surface we consider. For all values of A 7^ 0, 00, Xj there are two points on 
the Riemann surface, one in the upper sheet, and one in the lower sheet. 

Consider now Wj=i... 9 to be the unique basis of holomorphic abelian dif- 
ferentials satisfying § a ujj = 5^, and define the g x g period matrix 



(3.45) 



It can be proved that II is a symmetric matrix and its imaginary part is 
positive definite allowing the definition of an associated ^-function 



e(Q = ^ e 27ri ( n * nn+ " tc ) . 



(3.46) 



n6» 



The arguments of the ^-function are ( which is a vector in C 9 and the period 
matrix IT (which we consider fixed and therefore do not explicitly write as an 
argument). The sum is done over all n G Z 9 , that is all order g vectors with 
integer components. All vectors (e.g. n,() are taken to be column vectors 
(and therefore their transposes n*, are row vectors). Simple but important 
properties of the theta function are 



o(-0 = o(C) , 

and the (quasi)-periodicity: 

e (c + a 2 + nAO = e -H A ^H^A9(() , 



(3.47) 



(3.48) 



where A 1; A 2 G Z 9 , namely are vectors with integer components. To shorten 
some equations, it is also useful to define the 6 function with characteristics: 



HO = 



A 2 



(C) = exp <^ 27TZ 



-A^nAi + -A*C + -A^A 2 
8 1 2 1 4 1 



0( C + ^A 2 



(3.49) 

where again A 1; A 2 e Z s . We introduced the notation 9 for this function 
because, in the rest of the paper, Ai and A 2 are fixed vectors. In particular 
from now on we are going to consider that A* A 2 is an odd integer which is 

also described as saying that ^ is an odd characteristic. In such case 

2 

we have 

0(-C) = e™^6(Q = -9(C) , (3.50) 
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as can be derived from the definition of 9 and we used that, in our case, 
A^A 2 is odd. In particular this implies 



(9(0) = 



9 



1 * 1 A 

-A 2 + -nAi 

2 2 



o , 



(3.51) 



namely the vector a = |A 2 + |nA! is a zero of the theta function. The 
(quasi)-periodicity of the theta function implies that 



e 



Ai + 2ei 
A 2 + 2e 2 



(£) = e inA ^ 2 9 



Ax 
A 2 



(0, 



(3.52) 



for any e lj2 G Z s . Therefore it only makes sense to consider A 12 modulus 
two, namely its components being zero or one. 

The most important property of the theta functions that we need in this 
paper is Fay's trisecant identity: 



0(0 U + + / = 71234 9 (C+ / w) ((+ / w) +71324 ^ (C+ / w) (( 

V Jp2 J P A / J P2 J P A J P A 

(3.53) 



with 



lijkl 



9(a + J p y)9(a + J v y 



pi 



(3.54) 



0(a + J p »% + J p » 

In these formulas are points on the Riemann surface, and a is a non- 
singular zero of the Riemann theta function, i. e. the function is zero but not 
its gradient. In particular we are going to use a = ^A 2 + ^IIAi which is a 

zero as noticed before. Also notice that the contour integrals f Pb Uj define a 

to Jp a J 

vector which from now on, following standard convention will be abbreviated 
as: 

ppb pPb 

(3.55) 



rPb rPb 
/ "j^ 

JPa JPa 



The function 7 may be viewed as a generalization of the cross-ratio function 
on CP 1 to functions on Riemann surfaces. Some immediate properties of 
these function are: 



71233 = 71134 = 1, 72134 = 7l234> 7l214 = = 7i 232 . 



(3.56) 



One important use of the Fay's Trisecant formula is that it provides a direct 
way of obtaining directional derivatives of theta functions or of ratios of them. 



h2 ' 
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Taking the derivative with respect to p\ and then letting p 2 — > pi we get 



D Pl In 



0(0 



P4 
P3 



-D pi \n 



% + L 



Pi- 
PS, ' 



ppiv 
>P4 ' 



D Pl 6(a)e(a + J p P :) ( C + J£ ) 0<C + JJ 



P4 
Pi 



.(3.57) 



Here _D pi indicates a directional derivative defined as (summation over j 
implied): 

D pi F(C)=u^0^, (3-58) 

and should not be confused with a derivative with respect to pi that, if 
appears, we will denote as d pi . Also, the final expression is simplified using 
the identities (3.56). We can further derive with respect to p% and take 

D Pl 9(a)D P3 e (a) 9 (C + J£ ) # (C + J ? 



P4 - > P3 obtaining: 
D P3Pl In 0(C) = -D P3Pl ln0 I a + 



pi 



P3 
Pi 



rp3 
'pi 



02(0 

(3.59) 

This summarizes the basic properties we need. Much more is known about 
these functions as can be found in the references [21], 1 201. 



3.2 Solution to cosh-Gordon equation 



Eq.(3.59) shows that the second derivative of the logarithm of a theta func- 
tion contains the theta function. So solutions of 



dda = 2cosh2a = e 2a + e~ 2a 



(3.60) 



should naturally be sought as logs of theta functions. To eliminate the con- 
stant term in eq.(3.59) we subtract two such derivatives and get 

o(0 



D PlP3 In 



rp3t 
>Pl ' 



D pi e{a)D P3 e{a) 



(3.61) 



0^+2^)0(0 6>(c+j£-2j£Mc+j; 



rps- 

'Pi ' 



0(a + /;;)0(a + /-) i nc + C) 



^P3^ 

'Pi 



P3 * 

Pi 



2 (C) 
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To get back the same theta functions we need to exploit their periodicity and 
therefore require 



in 



to 



a 2 + nAi , 



(3.62) 



pi 



where A 2 , Ai are integer vectors and II is the period matrix in eq.(3.45). 
This gives 



6(( + 2 



in 



-2niA\(-iiTA\nAi 



0(C). 



pi 



We obtain 
n i g (0 

™ 0{C + L 



P3\ 
pi ' 



X 



e -*riA*C e -f A*IIAi 



D Pl d{a)D p ,d{a) 

32/ 



-i§ A\UA 1 



-e 2 



o\0 



P3^ 
Pi ' 



^7rA t 1 A 2p 27riA t 1 C e f A*nAi 



(3.63) 

(3.64) 



2 (O 



We should now choose pi, and the path of integration between them such 
that A[A 2 is odd so that e i7rA i A2 = -1. Then we take 



„2a 



-27riA t 1 C-^ L A t i nAi 



! (C) 



o\0 



and 



C = 2uj(p 1 )z + 2u(p 3 )z. 



(3.65) 



(3.66) 



The last choice results in d z ( = 2D P3 (, d( = 2D Pl (. The correct normaliza- 
tion for ( and a follows from the result 



D Pl 9(a)D P3 6(a) 
9 2 (0) 



if A|riAi 



(3.67) 



which is explained in the appendix. In any case it should be clear at this 
point that the overall normalization of a can always be adjusted so that 
eq.(3.64) becomes the cosh-Gordon eq. dda = 2 cosh 2a. The final and very 
important point is that the theta functions are generically complex but a 
should be real. Again, following [H] we impose a reality condition as follows. 
Suppose there is a g x g symmetric matrix T such that 



n = -rnr, c = -T(, 



(3.68) 
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Then, it is easy to prove, using the definition of the theta function that 0(0, 
0(C) are real whereas for example e i7rA ' c 0(C+ J£\j) = e^ A i^(C+|A 2 + ^nAi) 
is purely imaginary. This explains the minus sign in ( 3.65[ ) and proves that a 



is real. The matrix T is constructed [TT] from an involution of the Riemann 
surface that shuffles the basis of cycles (a^bi). To prove for example that 
6(0 is real we use 

6(0 = e - 2ni{nt ^ ntfln) (3.69) 
= ^ e - 2 ^(-™ t ^-5 n ' TnT ™) (3.70) 

= 0(C), (3.71) 

where in the last equation we redefine the summation variable n — > Tn and 
used T* = T, T 2 = 1. 

3.3 Solution to equations for ip, ip 

In the previous section we showed in detail how to use the properties of the 
theta function to solve the cosh-Gordon equation. Now we are going to do 
the same for the equations determining ip but in a more sketchy way. Notice 
that ip and ip are two linearly independent solutions of the same equations: 

(3.72) 
(3.73) 
(3.74) 
(3.75) 



dipt 


= --daipi 


+ \e a ip 2 


dip 2 


= e- a iPx + 


-daip 2 , 


dipi 


= ^Baip! - 


• e"> 2 , 


Bip 2 


= \e a ip Y - 


^daip 2 , 



which are the expanded version of eq.(2.42). As a first step we can define a 



function F = e a p- which satisfies 



8lnF= \e 2a ^= - e~ 2a F. (3.76) 
A F 
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By using the identities for the first derivatives of the theta functions and the 
value of e 2a already given one can readily see that 

2D P A«M« + Q __ 2 , a!c _ ¥A!nA , 0(OHC + Q ... 

where we introduced another special point in the Riemann surface that we 
call p4. The points p\$ are going to be taken as branching points, in particular 
for definiteness we take pi = and p 3 = oo. On the other hand p 4 is not a 
branching point and we take it to be on the upper sheet with p 4 = A, the 
spectral parameter. Going back to the equations for ^i j2 and using the same 
techniques we find that 

, r -U 6 ^ + Ipi + fpi ) 2zD M lng(ffi)+2z£> P1 ln0( O +ff. 1 )+2ir<gA«h;(0) 

^ *(C + J?) 

the constant C is determined to be 

2A,,0(a)0( a -Q , nA| 

C " «(£)«(») e ■ <3 ' 79> 

Again, we emphasize that the technique is to match the equation with the 
properties of the theta functions and choose the parameters appropriately. 
Another, linearly independent solution can be obtained by choosing a differ- 
ent point p4 that we call p 4 . However it has to be associated to the same 
value of A and therefore it can only be the same point but on the other 
(lower) sheet of the Riemann surface. Namely both p 4 and p 4 project on A. 
It should be noticed that, when pi = 0, namely one of the branching points, 
this implies 

/ Uj = - / Uj , (3.80) 

Jpi Jp! 

because the first integral is done on the upper sheet and the second one on 
the lower sheet where the function fi changes sign (and therefore Uj changes 
sign). We obtain 

A — . 1 . , ( f t V ( 3 - 81 ) 
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The (constant) normalization factor can be computed using the trisecant 
identity to give 



2D P3 9(0) 



-2D P3 9(a)e*-i — e 2 (j Pi ) ' 



Pi ' 



(3.82) 



(3.83) 



To finish this section we rewrite the solution using the function 9 to obtain 

PP4\ 



^2 



m 0(O Ho 



»(C) 



with 



\i = -2D p:j \n9( 



PA 



v = -2D P1 \n§( 



pi 



in 



(3.84) 
(3.85) 
(3.86) 



pi 



It is straight-forward to check directly that these functions satisfy equations 

(3.87) 



(3. 72), (3. 73), (3. 74), (3. 75). The only identities that are needed are 
- W Pl 9(a)D P3 9(a)e^ A * nAl = AD Pl 9 (0)D P3 (9(0) = #(0) 2 



and 

A = — 4e 1 1 J ' 



i J 4 
PI 



D K 9{a)9a 



P3 ' 



0(a + J p p ;)9(0) 



Pi 
pi 



-i 2 



*(J£)*(0) 



(3.* 



which are explained in the appendix. The last identity allows us to define 
(that is to appropriately choose the sign of the square root) 



-A = 2- 



rp4\ 
>Pl ' 



Then the final form for ipi 2 is simply: 

He + 1 



'^2 



0(0 



(3.89) 

(3.90) 
(3.91) 
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where the sign of the square root is chosen according to the previous equation 



and fx, v were defined in eq.(3.86). We can also compute (remembering that 
J Pi = — J Pi because pi = 0.) 



rP4 



9(0 

0(C - D i 

^ = Q( ,\ P 1 e* a e-»*- v ~ z . (3.93) 

At this point we can replace fa^ and fa^ in A and then in X. This allows 
us to compute the solution directly in Poincare coordinates as: 

*(2 JZ 



>P\ > 



/1Z + UZ I 



\mmm\\e 
\Hc-n\ 2 +\d(c-n\ 2 



(3.94) 



r p4^ 



etc - r p4 )6(( + p 4 ) - ere - r wc + r , 

I^C-OP + I^C-OI 2 
3.4 Shape of the Wilson loop 

The shape of the Wilson loop is determined by the intersection of the surface 



with the boundary. The boundary is located at Z = which, from eq.(3.94) 
and for finite z, z only happens if 

Z = 6(C) = or 9(C) = 0. (3.96) 

Later on we are going to deal just with the second case so we determine the 
shape of the Wilson loop by 

0(C) = 0. (3.97) 

This equation defines a curve in the world-sheet which in turn is mapped to 
a curve in the (X, Y) plane using the solution to the equations of motion 



(3.95) 



3.5 Computation of the Area 

The expectation value of the Wilson loop is determined by the area of the 
minimal surface we described. In conformal gauge the area is computed as: 

A = 2 J dX^dX^dadr = 2 J Adadr = A J e 2a dadr , (3.98) 
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where we used eqns.(2.6 ) to write the area in terms of the Lagrange multiplier 

A = 2e 2a , (3.99) 



A and then used that 



as follows from eqns.(2.21) and (2.37) using that |A| = 1. We could in 



principle replace a by its expression in eq.(3.65) and evaluate the integral 



numerically but such procedure fails because the area is divergent. The 
correct procedure is to identify the divergent piece analytically and then 
extract the finite piece in terms of a finite integral that can be easily evaluated 
numerically. In order to do so we use eqs. (3.65) and (3.59) (using a = 



to get: 



4{D piP3 \n9(0)-D plP3 In 0(C) 
4D piP3 In 0(0) -00 In 0(C). 



(3.100) 
(3.101) 



The first term in the expression for e 2a is a constant and the second one 
is a total derivative. The first integral is clearly finite but the second one 
contains the divergent piece that we need to regulate. In order to do that we 
observe that for these solutions 



Z=\9(C)\h(z,z) , 
where Z is one of the Poincare coordinates and 



h(z, z) 



P4\ 
Pi ' 



P4\ 
Pi ' 



|0(O)0(C)| \e 



fJ,Z-\-l/Z I 



m - o\ 2 + m - 1::)\ 2 



(3.102) 



(3.103) 



pi 



Furthermore using Stokes or Gauss theorem we find that for any well-behaved 
function F: 



dadrddF 



dadrV 2 F 



n ■ VFdi 



(3.104) 



where the contour integral is over the boundary of the Wilson loop (in the 
world-sheet), n is an outgoing normal vector and dl is the differential of arc 
length. The area is then: 

A = l6D PlP3 In 0(0) J dadr + j h-V\nh dl - j h-V\aZ di. (3.105) 
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The last integral is divergent and we concentrate now on extracting the lead- 
ing divergence. The correct AdS/CFT prescription is to cut the surface at 
Z = e and write the area as 

A = ^ + A f , (3.106) 

where L should be the length of the Wilson loop and Af is the finite part 
which is identified with the expectation value of the Wilson loop through: 

(W) = e~£ A f , (3.107) 

where here A is the 't Hooft coupling of the gauge theory (not to be confused 
with the spectral parameter). This prescription is equivalent to subtracting 
the area A = — of a string ending on the contour of length L and stretching 
along Z from the boundary to the horizon. To see that the coefficient of the 
divergence is indeed the length, let us compute 

A dw =-<£ -h ■ VZdi = - I \VZ\di , (3.108) 

JZ=€ % e JZ=e 

where we used that the normal is precisely in the opposite direction of VZ 
because the contour is a curve of constant Z = e and Z increases toward the 
inside. On the other hand the length in the boundary is given by 



|t.VX| 2 + \t.VY\ 2 d£ , (3.109) 
where t is a unit vector tangent to the contour. We can move forward if we 



write the equation of motion for X as derived form the action (2.4): 

2VX-VZ = ZV 2 X, (3.110) 

which, when Z — y 0, becomes VX ■ VZ = namely VX is perpendicular to 
the normal and therefore parallel to the tangent t. The same is true for Y 
so we find 

L = j VIVXI 2 + \VY\ 2 d£ + 0(e 2 ). (3.111) 
Finally the equation of motion for Z is 

(VZ) 2 - ZV 2 Z = (VX) 2 + (VF) 2 , (3.112) 
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which for Z ->• implies that a/|VX| 2 + |VF| 2 = |VZ|. Therefore the 
length of the Wilson loop is given by 

l = j\ vz \ M -\j^Y\ M ' ^ 3 ' 113 ) 

and the divergent piece of the area is indeed A^ v = j. There is a finite 
part remaining: 

A = ^ + A { , (3.114) 



A { = 16D PlP3 In 9(0) J dadr + jn ■ Vln/i di + ^ j> 



V 2 Z 
VZ 



The integrals are performed on the world-sheet parameterized by a, r. The 
first integral is proportional to the area of the world-sheet. The last two 
integrals are done over the world-sheet boundary. The final expression can 
be simplified by rewriting Z = \9(()\h(z, z) and using that 9(Q vanishes on 
the boundary where the contour integral is performed. It is then easy to 
check that h(z, z) cancels and the final formula for the area is: 

A = ^ + A { , (3.115) 

A i = 16D piP3 In 0(0) J dadr di (3.116) 

= 16D P1P3 In 61(0) ! dadr - 2 I di (3.117) 

= 8D pip:i In 0(0) l{adr - rda) - 2 <f Dp ^ 9 ^ di. (3.118) 
J J \D P1 9(C)\ 

where in the last step we used the well-known formula for the area of the 
region encircled by a given curve. The renormalized area is then expressed 
as a finite one dimensional contour integral over the boundary of the world- 
sheet. Perhaps it would also be useful to clarify that \V9(Q \ denotes the norm 
of a real 2-vector whereas \D Pl 9(() | is the modulus of a complex number. The 
final expression for the area is quite interesting because it does not depend on 
the spectral parameter A. Therefore the shape of both, the Wilson loop and 
the dual surface depend on the parameter A but the area Af does not. In this 
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A 

T 




C 

Figure 1: The boundary is determined by the contour Z = 0. However the 
area is computed by integrating up to a contour Z = e — > and then the 
leading divergence - is subtracted. Here L is the length of the contour in 
the boundary (not in this (cr,r) plane). 



way we explicitly find a one parameter family of deformations that preserve 
the area. It is not obvious at first that the area Af should be independent 
of the spectral parameter because, although the definition (3.98) does not 
contain A, the regularized area does since, as we said, L depends on A. It so 
happens that the finite part Af does not. The situation is similar to scale 
transformations that modify L but not Af. 



4 Example with g = 3 

We are going to consider an example to illustrate the shape of the Wilson 
loops that are obtained in this way. The main purpose is to show that we 
find closed Wilson loops whose dual surface is known analytically. 
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Consider the function 



/i = iy/-i(\ + 1 - i)y/-i(\ + 1 + i)y-«(* - - ^)v / 2^Av / xy A + i 

(4.119) 

where the square root is taken to have a cut in the negative real axis. So 
defined, the function \i has cuts in the complex plane as illustrated in figj2] 
but is smooth in a double cover of the plane which defines a hyperelliptic 
Riemann surface of genus g — 3. The cycles a*, 6, are taken as in the figure. 
The Riemann surface has an involution A — > — j meaning that knowing 
the cuts for |A| > 1 we can reconstruct the cuts inside the unit circle. This 
involution is important to construct the matrix T that fixes the correct reality 
conditions (see [11J). A basis for the holomorphic abelian differentials is given 
by 

\k-l 

v k = dX, A; = 1 ... 3. (4.120) 

fi 

If we compute 

Cij = j> Uj, Cij = j> Vj , (4.121) 

J Obi *J °i 

then a normalized basis of holomorphic abelian differentials is 

uh = {C^) jt , (4.122) 

and the period matrix is 

0.5 + 0.64972i 0.14972i -0.5 
U = CC- 1 =\ 0.14972i -0.5 + 0.64972i 0.5 |. (4.123) 

-0.5 0.5 0.639631 

The Jacobi map (with base at 0) is defined as 

0(A), = f ^ = f" — (C- 1 ) dX. (4.124) 
Jo Jo f 1 

The function 6(4>(X)) has three zeros: A = oo, i=i,-l + i. To prove this we 
can take, on the upper sheet, a path from A = to each of the zeros of fi. 
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Coming back along the lower sheet defines a closed path C\ equivalent to: 



A = 


i 

2 


— y 


C A 

^ A 


= b?. — a? 


— >■ 





A = 


- - -i 

2 2 L 


->• 


C x 


= a 2 + b 2 


-)■ 




A = 


2 ^ 2 l 


-)■ 


Cx 


= b 2 


— > 




A = 


2 


— >■ 


Cx 


= a 3 — a 2 


-> 




A = 




-> 


Cx 


= -a 2 + a 3 + 6i + 6 3 


-> 




A = 


-1 - i 


-)■ 


Cx 


= -a x + a 2 + a 3 + 61 - 63 


->■ 




A = 


OO 


-)■ 


Cx 


= ai — a 2 + a 3 + 63 


->■ 






(0 
(0 



1 
0-10 
1 
1 
1 
000 
000 

-1 1 

1 1 
-1 1 

1 -1 
-1 1 1 

1 

1 -1 1 

(4.125) 

The last column gives a theta function with characteristic determined by 



(0 



(0 
(0 

(0 
(0- 



£1 
£2 



where £i j2 are given by Cx = 2^i=i( £ iA + e 2i a i)- The point is that 

this theta function vanishes if and only if #(0(A)) vanishes because 0(A) = 
|(e 2 + Ilex) (since the integrals on the upper and lower sheet are equal to 
half the total integral). If the characteristic is odd (namely e\e 2 is odd) then 
by symmetry the theta function is zero at the origin and therefore the theta 
function without characteristic is zero at the corresponding point. 

The vector of Riemann constants [21 J is computed from the sum of the 
characteristic of the zeros which is (mod 2) 



K 



I 1 
1 1 








(4.126) 

and has the property that for any Ai 2 on the Riemann surface we have 







1 1 
1 1 








(0(A 1 ) + 0(A 2 )) = O. 



(4.127) 



Moreover, this completely defines the set of zeros of the theta functions we 
are considering [2T] . 

To write the solution we choose the points pi — and p% = 00 which 
works well since, as seen above a path between and 00 defines and odd 
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a, 



b, 



a 3 



Figure 2: Genus three Riemann surface. A basis of fundamental cycles is 
depicted with solid paths being in the upper sheet and dotted lines in the 
lower sheet. The circles show their intersection points. 
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1 




characteristic that can be taken to be Ai = 





, A 2 = 


1 


. We there 




1 




1 





fore define 



9(C) = 9 



1 

1 1 1 



(0- 



We can now write 

C = 2(u(oo)z + u{0)z) 



-0.4903a - 0.19069r 
-0.4903(T + 0.19069r 
0.59321r 



(4.128) 



(4.129) 



where we defined z = a + it. The zeros of the function 9 in the complex 
plane z determine the boundary and therefore the shape of the Wilson loops. 
We plot the contours where 9 becomes zero in figji 

The (quasi)-periodicity of the theta function is evident from the figure but 
also the existence of closed curves which in turn will give rise to closed Wilson 
loops. Choosing the contour displayed in fig|4] determines a Wilson loop up 
to the spectral parameter A that is still arbitrary. Taking into account that 
|A| = 1 we choose as examples 



Ai 



A 2 



1+j 
v2 ' 



(4.130) 



in the formulas for the surface 



Therefore the point p 4 = i (or p 4 = — ^] 
(understood to be in the upper sheet). The shape of the Wilson loop can 
be immediately obtained by mapping the contour into the X, Y plane. The 
result in displayed in fig |5] where we see a non-trivial closed Wilson loop. The 
minimal surfaces ending on these contours are displayed in fig|6j Finally we 
can compute the length of the Wilson loop and the finite part of the area 
using eq.(3.117). The result is 



Lt = 13.901, 
A f = -6.598 



L 2 = 6.449 
for both. 



(4.131) 
(4.132) 



The length can be changed by a scale transformation but the finite part 
Af is scale invariant and independent of A, the spectral parameter. For 
comparison, for a circle of radius R the corresponding values are: 



L (circle) = 2?r d A 



(circle) 
f 



-2vr. 



(4.133) 
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-2-10 1 2 3 

IT 

Figure 3: Zeros of the function 9. We can see the (quasi)-periodicity of the 
theta function as well as closed contours that map into closed Wilson loops 
in the boundary. 
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2.0 - 



1.8 - 



1.6 - 



1.4 - 



1.2 - 




0.0 0.2 0.4 0.6 0.8 1.0 

(T 

Figure 4: Particular contour (taken from flgj3]) in the z plane chosen to 
compute the minimal area surface. 
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(a) A = i 



Figure 5: Shape of the Wilson loops that we use as an example. It is obtained 
by mapping the contour in fig|4] into the X, Y plane for two values of the 
spectral parameter A. 



Although in the end the results are numerical, we emphasize that the shape 
of the surface is known analytically. Also a relative simple expression was 
found for the expectation value (after subtracting the infinities) in terms 
of a one-dimensional finite integral. At the end, the integral was evaluated 
numerically to get the value for the area. Notice also that the area is smaller 
than the one for the circle. For a given length L the area is not bounded 
from below since we can take a contour made out of two parallel lines of 
length L/2 and separated by a distance 5 — > in which case the area goes 
to minus infinity as Af ~ — ^ [2j. On the other hand, for fixed L, the circle 
is expected to be an upper bound as shown in [22]. Our result agrees with 
that bound. 

Having described a particular example in detail we want to elaborate 
further on the properties of these Wilson loops. We know that th e zero s 
of 8(() determine the boundary of the Wilson loop, but, from eqs.(4.127), 
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(a)A = i ( b ) A = -7? 

Figure 6: Minimal area surfaces ending on the contours illustrated in figj5} 
We emphasize that the surfaces are known analytically. 



(4.128) we know that all zeros are given by 



" " 




" 1 " 







+ in 

2 


1 


+ 0(Ai) + 0(A 2 


1 




1 





(4.134) 



for arbitrary points A12 on the Riemann surface. The vectors 





1 



and 



9 



represent the difference (mod 2) between the characteristics of 9(0 
(() and k = 



1 

1 1 1 



1 1 
1 1 



on the Riemann surface because we a 

( = 2 (UooZ + u z) , 

which satisfies 

0-10 

<, = t = \ -1 

00-1 



. We cannot take random points Ai^ 
so have 



(4.135) 



(4.136) 
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With some work it can be seen that we can take 



2 m 




" 1 +2 mi " 


2n 2 


+ in 


1 + 2m 2 


l + 2n 3 


2 


_ 1 + 2 m 3 _ 



-0(Ai) + 0(-i- 



(4.137) 



where m,2,3 and mi )2 ,3 are integers which can be absorbed in the definition of 
the path used to compute the function 0. It therefore follows that for genus 
3 we can map the Wilson loop into a curve inside the Riemann surface. 
Namely for each point in the contour displayed in fig|4] there is a point Ai in 
the Riemann surface such that: 



2(u(oo)z + u)(0)z) 



1 
2 



2n x 




" l + 2m 1 ' 


2n 2 


+ in 


1 + 2m 2 


l + 2n 3 


2 


_ 1 + 2 m 3 _ 



-0(Ai)+0 



1 



(4.138) 

The set of such points describes a curve inside the Riemann surface which 
is depicted in figj7| The statement is the following: for each point Ai in the 
curve and for a given choice of path used to define the function cf> there is a 



set of integers ^1,2,3 and mi )2 ,3 such that eq.(4.138) can be solved for z. These 
values of z lie in the closed curves depicted in fig|3] which are the zeros of 
the function 9. Equivalently we can set all integers ^4,2,3 = and 777,1,2,3 = 
but then we have to choose the path used to define appropriately so that 



eq.(4.138) has a solution. Furthermore, with an appropriate choice we can 
map the curve in figjT] to the curve in figjl] which was the one used in our 
examples. For higher genus, there should be a set of curves in the Riemann 
surface for each closed Wilson loop. In this paper we do not explore this 
issue further but we believe it should be an interesting subject to pursue. 

As a final comment, it should be noted that these Wilson loops are not 
BPS since, for Euclidean Wilson loops with constant scalar, the only BPS 
ones are straight lines 



5 Conclusions 

In this paper we discussed minimal area surfaces in AdS^ space which are 
dual to Wilson loops in = 4 SYM. We essentially follow the results of the 
paper [UJ where such solutions were found but provide a different derivation 
of the solutions and also a formula for the finite part of the area in accordance 
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Figure 7: For genus three, each Wilson loop computable in terms of theta 
functions maps into a closed curve inside the corresponding Riemann surface. 
In this figure we display the closed curve corresponding to the Wilson loop 
of fig|5j We notice it encircles two cuts. 

with the one used in the AdS/CFT correspondence. Finally we make the ob- 
servation that closed Wilson loops appear from these solutions in which case 
the world-sheet has the topology of a disk and the area is expressed as a 
(finite) contour integral over the world-sheet boundary. In this way we con- 
struct an infinite parameter family of new examples of Wilson loops whose 
dual surface is analytically known. It should be noticed that, to our knowl- 
edge, only the circular Wilson loop and the lens shaped loop were previously 
known for closed Euclidean Wilson loop (with constant scalar). Furthermore, 
the result gives, for each individual Wilson loop a one parameter family of 
deformation (given by the spectral parameter A) such that the area remains 
the same. Finally, for genus 3 we pointed out an interesting map between 
Wilson loops computable in this way and curves inside the Riemann surface. 
We hope that these new solutions will give rise to a better understanding 
of Wilson loops in the context of the AdS / CFT correspondence and also in 
general. It is evident that an important integrable structure lies behind them 
that should be explored in detail. It would be interesting if these ideas allow 
us to reconstruct the shape of the Wilson loop from the field theory using 
some sort of coherent states formalism as in |24J. 
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7 Appendix 



In this appendix we derive a useful identity for the theta functions. Instead 
of doing a general derivation we show how it works in the example we are 
dealing with in the main text and then the generalization should be clear. 
Consider the function 



tA* r 
i j? 



h(X) 



0{a + f: 



(7.139) 



pi ■ 



defined on the Riemann surface. Also c is a constant that we choose later. By 



looking at the results (4.125) one can see that the numerator has zeros at A = 
0, i), whereas the denominator vanishes for A = oo, i), 

It follows that h(X) has a zero at A = and a pole at A = oo. To see the 
behavior near A = we use that (remembering that pi = and eq.(4.122)) 



J Uj ~ J Tj^ c i 3 l = -2zVACy 1 , A -> 0. 



Thus, 



If we choose 



h(X) ~ -4cA 



D pi 6(a) 



9(0) 



c = — 



1 / 9(0) 



4 \D pi 9(a) 



(7.140) 



(7.141) 



(7.142) 



then h(X) ~ A, that is, it has a simple zero at A = 0. Similarly one can check 
that it has a simple pole at A = oo. Finally it can be seen to have the right 
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periodicity properties to be well defined on the Riemann surface. The only 
such function is h(X) — X so we conclude 



A _ 1/ 0(0) V e'^e(a + Q 

x -^{D^m) [ ^ 1 ■ (7143) 



as used in the main text. Furthermore, by expanding around A = oo we find 
that h(X) = X if we choose the constant c to be 

c = _ 4e - A i nA 1 ( 2 . (7.144) 



9(0) 

Equating the two expressions for c we find 



D Pz 9(a)D Pl 9(a) \ 1_ i7rAi n Al 



2 

' e - i7rA i nAl , (7.145) 

as we also used in the main text. Although we derived the result for our 
particular case it is clearly valid in general (see for example [21] for the case 
of a generic hyperelliptic Riemann surface). Taking the square root of the 
last equation we find 

Dp 3 9(a)D Pl 9(a) _ 1 ., A;nAl , , 

The sign cannot be determined by this reasoning but it is easily found to 
be minus by a simple numerical computation. The result was used in the 
main text to find the correct normalization of ( so that a is a solution to the 
cosh-Gordon equation. 
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